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COOLING OP GAS TURBINES 
TV - CALCULATED TEMPERATURE DISTRIBUTION IN 
THE TRAILING PART OP A TURBINE BLADE 
USING DIRECT LIQUID COOLING 
By W. Byron Brown and William R. Monroe 


SUMMARY 

A theoretical analysis was made to determine the temperature 
distribution in the trailing portion of a liquid-cooled turbine blade 
between the coolant passage and the trailing edge and to study the 
Influence of various design and operating variables on the hot-Bpot 
temperature at the trailing edge. The trailing portion of a typical 
turbine blade was approximated by equivalent rectangular and wedge- 
shaped sections to facilitate analysis, and representative values 
for the boundary temperatures and surface heat-transfer coefficients 
were assumed or computed according to standard methods. 

Three-dimensional temperature distributions were obtained at 
effective gas temperatures of 2000°, 3000°, and 5000° P for a 4-inch 
prismatic section, cooled over the entire length and having a thermal 
conductivity of 15 Btu per hour per foot per r, and at 2000° P for 
similar sections having thermal conductivities of 120 and 210 Btu 
per hour per foot per °P. In addition, comparison of one-dimensional 
temperature distributions through the outer portion of the blade was 
made for two shapes, equivalent rectangular and wedge plane sections, 
at an effective gas temperature of 2000° P« The influence of thermal 
conductivity on one-dimensional temperature distribution in both 
plane sections was evaluated. With the equivalent rectangular solu- 
tion, an investigation was made to determine the effect on trailing- 
edge hot-spot temperature of variations in the trail? ng-section 
width, thickness, thermal conductivity, and coolant-passage heat- 
transfer coefficient. 

The three-dimensional and one -dimensional solutions of the 
equivalent rectangular trailing portion give almost Identical temper- 
ature distribution in the blade section away from the influence of 
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rim cooling. For the "blades of low thermal conductivity and wide 
trailing sections the hot-spot temperature closely approaches the 
effective gas temperature, and infinite increase in the surface heat- 
transfer coefficient at the coolant passage relative to that obtained 
■with water will have small effect on the trailing-edge temperature. 
Further reduction of trailing-edge hot-spot temperature can be 
obtained by the following three methods given in order of decreasing 
effect; by use of materials having high thermal conductivity, by 
altering the blade shape to provide a relatively short trailing 
portion, by locating additional cooling passages a short distance 
from the trailing edge. 


INTRODUCTION 

As part of the general program to increase the permissible 
operating temperature of gas-turbine cycles and to improve the life 
of critical parts, further analysis of the liquid-cooled turbine 
blade is being conducted by the MCA to determine the available 
cooling effect at the anticipated hot spot under present conditions 
of design and to evaluate various possible methods of Improving 
cooling of the blades. 

Part I of the current Beries on "Cooling of Gas Turbines" (ref- 
erence l) presents calculations on the application of air cooling 
fins to the turbine rotor with the object of reducing. the rim temper- 
ature to permit better cooling of solid blades. It was found that 
the root of the blade could be cooled about 800° F below the effec- 
tive gas temperature, but relatively little gain was obtained at the 
critical section some distance from, the root. In part II (refer- 
ence 2) it was shewn that blade life could be extended by rim 
cooling, and that the effective gas temperature could be slightly 
increased . 

Part HI (reference 3) discusses the circulation of liquids 
through hollow passages in the blade and presents a one-dimensional 
analysis of temperature distribution from the blade tip to the rotor 
axis. The effects of variations in blade length, length of cooling 
passage, type of liquid coolant, and rate of coolant flew were 
studied. The temperatures calculated (one-fifth the gas temperature 
for water cooling) were those appropriate to points near the cooling 
passages. 

The purpose of the present detailed study is to investigate the 
temperature distribution encountered in the relatively long, thin 
trailing section "between the coolant passage and the trailing edge 
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of the liquid-cooled "blade. This section, particularly the trailing 
edge, seems most likely to "be critical. The "blade design and 
cooling conditions used are the same as used in reference 3 except 
that only the "blade section bounded by the rim, the tip, the wall 
of the coolant passage, and the trailing edge is considered. The 
analysis is carried out in several parts, with assumed shapes or 
sections that approximate the trailing portion of a typical blade 
as follows: a one -dimensional analysis of a rectangular section 

through the outer portion away from the root, a similar one- 
dimensional analysis of a wedge-shaped section, and a three- 
dimensional analysis of a right prism for the entire span of the 
blade. 

The following sections of this report present the development 
of the various solutions for temperature distribution in terms of 
the dimensions, physical properties, boundary temperatures, and 
surface conditions that prevail. In addition, a number of figures 
are presented that illustrate the various effects of thermal 
conductivity on three-dimensional temperature distribution for 
several values of the effective gas temperature and the combined 
influence of design dimensions, thermal conductivity, and some 
surface conditions on cooling of the trailing edge. 


THEORETICAL ANALYSIS 
Symbols 

The following symbols, arranged in alphabetical order for the 
convenience of the reader, were used in the calculations: 

a 

b 

k 

l 

*i 


t 



■parameter eq.ua! to > (ft“^) 

distance from trailing edge to wall of coolant passage, (ft) 

thermal conductivity of turbine metal, Btu/(hr) (ft) (°F) 

distance from blade tip to rim, (ft) 

heat-transfer coefficient from hot gas to metal, 

Btu/(hr)(sq. ft) (°F) 

heat- transfer coefficient from metal to coolant, 

Btu/ (hr) (sq. ft) (°P) 

blade thickness, (ft) 
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t mean blade thickness, (ft) 

ti blade thickness at trailing edge, (ft) 

tg blade thickness at wall of coolant passage, (ft) 

T temperature of metal, °F 
T e effective temperature of hot gas, °F 

average temperature of liquid coolant, °F 
T r temperature at rim, °F 

x distance from blade tip to blade element, (ft) 
y distance from trailing edge to blade element, (ft) 
z distance from median plane of section to blade element, (ft) 

2a angle between sloping sides of wedge 
Q T 0 - T, °F 

A linear dimension accented by a prime denotes that the distance 
has been extended by tj/2 (for the wedge) or t/2 (for the rectangle) . 


Discussion of Simplifying Assumptions 

The turbine wheel with cooling passages is shown in figure 1 
and a cross section of the blade in figure 2. The approximate blade 
shapes (right prism, wedge, and rectangle) used in this analysis, which 
can be used to approximate any tapered profile, are given in fig- 
ures 3, 4, and 5. 

The equations for the temperature distribution in the blade 
section between the coolant passage and the trailing edge were 
obtained by equating the heat entering and the heat leaving an 
element. In order to simplify the boundary conditions, the cross 
section of the blade under consideration was approximated by geo- 
metric figures having straight sides and areas equal to that of the 
blade cross section. For the first development, a rectangular shape 
(fig. 3) was used, in which the width was made equal to the distance 
from the trailing edge to the coolant-passage wall and the thickness 
eq.ual to the mean thickness of the blade over this width. For a 
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closer approximation to the true shape, a wedge shape (fig. 4) was 
used in which the width remained the same as in the rectangle and 
the thickness tapered from, the blade thickness at the coolant passage 
to the blade thickness at the trailing edge. For both shapes, the 
approximate cooling surface was less than that encountered in the 
actual blade because the circular coolant-passage wall was approximated 
by a chord rather than by the true arc. 

The analysis was performed in the following parts': (a) a one- 

dimensional analysis of a rectangular section through the outer 
portion of the blade, away from the blade root; (b) a one -dimensional 
analysis of a wedge-shaped section through the outer portion of the 
blade; and (c) a three-dimensional analysis of a right prism for a 
4-inch blade length. 

The following simplifying assumptions were made for these parts; 

(1) In parts (a) and (b), the heat flow to the rim. was negligible 
at a section chosen away from the blade root because the length l 
was large compared to the width b (fig. 5) . The validity of this 
assumption is demonstrated by the results obtained in the three- 
dimensional analysis* 

(2) The heat gained from the trailing edge can be accounted for 
by assuming the width to be extended by a distance equal to one-half 
idle blade thickness at the trailing edge t^ (fig. 6 ) . The edge 
OQ at temperature T gained same heat. The extended surfaces 00* 
and QQ* were at nearly the same temperature T and no heat entered 
the end 0 'Q ' , therefore these surfaces gained practically the same 
amount of heat as the actual exposed end. In part (c) especially, 
the boundary condition was much simplified. This assumption remained 
valid as long as t/2 was small compared with b. This validity is 
demonstrated for part (a) . In like manner, the blade length was 
extended a distance eq.ua 1 to t/2. The length and the width extended 
by the distance t /2 are denoted by b' and l ' , respectively. 

This assumption was used only in parts (b) and (c; . 

(3) The temperature of the blade at the rim was constant at rim 
temperature. If this constancy is not maintained, it is shown that 
variations from it have little effect on the results except very 
close to the rim. This assumption was used only in part (c). 

(4) The variation of k was negligible and T e , T} , < 3 . 3 . , and q 0 
are constant over the blade and surface passages. Because these 
blades were assumed to be cooled over the entire length and those of 
reference 3 cooled only to within one-sixteenth inch of the tip, the 
comparable curves differed at the tip. 
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The analysis of parts (a), (b), and (c) are given In the follow- 
ing sections. 


One-Dimensional Analysis of Rectangular Section through. 

Outer Portion of Blade, Away from. Blade Root 

The derivation of the formula for part (a) as shown by figure 3 
is obtained from: 

Heat entering element from right end = kt de/dy 
Heat entering element from, sides = 2q.^ dy & 

Heat leaving element from left end = kt (de/dy + d^/dy 2 dy) 
Therefore, 

kt de/dy + 2q. ± dy 0 = kt (de/dy + d 2 e/dy 2 dy) (l) 


or 


where 



( 2 ) 

(3) 


The boundary conditions are as follows : 
when y = b, 

- lot (T - T*) 

If T^ - 0 is substituted for T, 

tJ ‘ 


kdg 

dy 


<i 0 (T e - t i - e) 


when y=0, 




A solution of equation (2) is 


(4) 

(5) 


e = A cosh a (y + e ) 


( 6 ) 
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where A and c axe Integration constants. The quantities e 
and A are evaluated to satisfy the boundary conditions at the 
blade edge and in the cooling passage, respectively. 

Eram. equations (5) and (6), when y = 0, 

kA a sinh a £ = q.j A cosh a e (7) 

therefore from equations (3) and (7) 



when t is not too large. This difference between e and t/2 
varies from 0.4 to 0.5 percent for the values of a considered and 
has an inappreciable effect 0.C05 foot or more inside the trailing 
edge. 


When y = b, then from equations (4) and (6) 

*o 


Aa sinh a (b + c ) = — [T e - Tj - A cosh a (b + 6 ) ] 

k 


or 


A = 


IT (T e - ? z ) 


a sinh a (b + e ) + ~ cosh a (b + e ) 

k 


( 9 ) 


( 10 ) 


Therefore, 


e = 


q__ 

-£• (T© - T Z ) cosh a (y + 6 ) 

a sinh a (b + c ) • + cosh a (b + e ) 

ic 


( 11 ) 


One-Dimensional Analysis of Wedge-Shaped Section 
through Outer Portion of Blade 

It has been demonstrated (reference 4, equation (53)) that a 
solution for the temperature distribution along the y-aiis in a 
wedge-shaped segment (fig. 4) may be expressed 

6 = A'J 0 (in) +BiE 0 (ip) (12) 

where A' and B are arbitrary constants and J D and Hq are 
Bessel's functions of two kinds, zero order. 
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The variable [j. is defined 


H 2 = 4C 2 


/’ + *1 


f 1 - tan a \ 
2 tan a ) 


where 


and 


c 2 = 


a = tan" 




k Bin a 


(t 2 - tA 

V 2b ) 


The boundary conditions are, -when y 1 = 0 

dO 


dy 1 


= 0 


when y ' = b ' 
or when |a - y.j 
and when p. = |a 2 
when y' = 0 
when y ' = b 1 


k ■— - Oo (To - - 0) 


M= 0 

d\i 


* < T e -*i-e> 

m - 2C 


- tan a) 


tan a 


^ - »f'**r(VTETr; 


(13) 

(14) 

(15) 

(16) 
(17) 


Equation (13) can bo differentiated and put in the form 


d0 _ _ji_ _d0_ 
da 2C 2 dy' 


de 

&\±* 


(18) 

(19) 

giving 

(20) 


Therefore when 


then 


dG_ 

dy' 

d£) 

d(-L 


= 0 


= 0 


By virtue of the properties of Bessel's functions, 

dJ Q (r)/dy ■ -J x (y) 


( 21 ) 
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and 


therefore 


and when 


then 


d6> 

d(i 


(y)/dy = -&!_ (y) 

= -A'iJj; (iji) +BE X (in) 


( 22 ) 


dn 


B = 


ia'Jx (im) 

H i Ufd) 


(23) 


In. order that equation (17) may he satisfied 
2 

h C”A'i J X (in 2 ) + BH-^iug)] = ^>[(T e -T 2 ) - A'J 0 (ip 2 ) - Bi Ho(in 2 ) ] 

(24) 

In standard form equations (23) and (24) can he written, respectively 
[«! (inx)] A« - IE ± (in x )] B = 0 


*r 0 (in 2 ) - Ui(in 2 ) 


r2C2fc 1 


A ’ + [5Tnr E l (i(1 2) + iBo(^ 2 )| B = T e -T 2 


<*0^2 


T3i© values of A 1 and B can easily he found and substituted into 
equation (12)^ giving Q in the form 

( 3joJ- l 2 

r (T e ~T 2 ) [HiUn^Jodn) + U 1 (i|.i 1 )iB 0 (in)] 

6 (25) 

[iJ 1 (iH 1 )H 1 (in 2 )3 - C3i(im)iJi(i|i 2 )] + 


■where X is 


CH^mJJodug) + iJxdHiJiBod^)] 


Three-Dimensional Analysis of Bight Prism for 
4 -Inch Blade Length 

The derivation of the formula for part (c) as shown in fig- 
ure 5 is obtained from 
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Heat entering element from top = -k dydz 


St 

51 


Heat entering element from right and = — k dxdz gjp 

ST 

Heat entering element from front = -k dxdy 

^St S 2 t 

Beat leaving element from hot ton = -k dydz 1 3 — + —sr dx 

\ ox Sx° 

/ 3t S^T \ 

Heat leaving element from left end = -k dxdz j ^ + g-g- dy I 

/St S 2 T \ 

Heat leaving element from hack = -k dxdy + g— dzj 


When the heat entering is equated to the heat leaving 

Sr' . . St , , . St 
31 - k dzdz 37 “ k dzd ^ 37 


, . , St ' . . St . , . St . , . /St S2t . \ 

■h dydz 3 ^ -k dxdz 33 : -k dxdy 33 - » -k dydz i gj + g-g- dx 1 


— k dxdz 


/St S 2 t . \ . . _ /St S 2 t \ f * 

[ 5 ? * aj2 ay ) - fc aiay var + az j <26) 


or 


0 _ SfT S2 t 

" Sx 2 Sy 2 Sz 2 


(27) 


The boundary conditions can he more simply handled if the gas 
temperature is taken as a reference temperature rather than the 
metal temperature; that is, if the substitution is made 


or 


Q = T e - T 


t = t 0 - e 


Equation (27) then becomes 


o - d 2 8 S 2 8 S 2 9 

ax' 2 Sy' 2 Sz 2 


(28) 

(29) 

(30) 


The origin of the coordinates chosen is indicated in figure 5. 
The axis z « 0 is located in the median plane of the right prism; 
from considerations of symmetry there can he no heat flow across the 
median plane. 
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The 'boundary conditions that must "be satisfied are 
when x‘ = 0 



ajar 

ii 

o 

(31) 

when y* = 0 

o 

II 

© 

row 

(32) 

when z = 0 


(53) 

when x* = 1' 

CD 

II 

►3 

CD 

1 

(34) 

when y* = b' 

n 

6 ° 

o* 

i 

H3 

1 

ct) 

(35) 

when z = 

• — - ■ 



kde 

5£“ = 0 

(36) 

A convenient form for the solution of equation (30) 

Is 

Ct CO 

e 

i i 



X (Arm COS I^X 

cosli 1 cos Pmz +Bmr> cosh. RmrPc 1 cos Apy 1 cos ! 



(57) 

where A, B, M, If, B, K } and P are arbitrary constants. The rela- 
tion between the constants used as coefficients for the variables Is 


f'W 2 = Pm 2 + N n 2 

(33) 


®nin 2 = ^m 2 + \n 2 

(39) 


tion of equation the"’ 'boundary condiiions 4 expreBsed 7 in^eaua- 

tlone (31) to (33) are satisfied* 


selec- 
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In appendix A, values of Amn and Bmn aro derived. These 
values are 

P m t 


(-l)n-l sln 8qo (T 0 - Tj) 


Affln it(2n-ij t 0TJ t Tdifem sinh Mmnb* + q 0 cosh Mmnb 1 

rjji n sin firm n 

1 + -T~ 

2P m | 


Bmn - 


sin \ n b» 4 (T e - T r ) sin F m jT 


An * 1 


1 + 


gpgil P„ | sin 2P m | 00611 1' 

***& * 1 + 

2P m f 


Values of n from 1 to 14 were Inserted and the temperatures cal- 
culated from equation (37). 

In appondix B it Is shown that sufficiently accurate results 
(within 1°) can he obtained by using only the first value of m; 
that is, m = 1 and dropping the terms involving cob Pg z, 
cos P~ z, etc. 


APPLICATION CO? ANALYSIS 

The following assumptions were made in applying the previous 
results to specific numerical calculations: 

1. The gas flow was 55 pounds per second. The heat-transfer 
coefficient corresponding to this gas flow qi was found to be 
222 Btu per hour per square foot per °3P. 

2. The liquid coolant, water, had an average temperature of 
200° P and a flow rate of 6.42 pounds per second for the entire rotor. 
The heat-transfer coefficient corresponding to this coolant flow q,., 

was calculated (reference 5) as 2370 Btu per hour per square foot per °P. 

3. Average thermal conductivities k of 15, 40, 60, 80, 100, 

120, and 210 Btu per hour per square foot per °P were used. 

4. The turbine had 55 blades. 

5. The blade was liquid cooled over its entire length of 4 inches 
by two cooling passages 0.25 inch in diameter._ The blade chord was 
1.188 inches long. 
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6. Some numerical values assumed were 
b = 0.0500 foot 

b* = 0.0550 foot (rectangular section) 
b* = 0.0515 foot (wedge-shaped section) 

Z = 0.3333 foot 
l* = 0.3383 foot 

t = 0.010 foot (rectangular section) 

t-j_ = 0.003 foot (wedge-shaped section) 

tg = 0.021 foot (wedge-shaped section) 

T r - 284° P (all values of thermal conductivity) when 
T 0 = 2000° P 

T r = 331° P (all values of thermal conductivity) when 
T e = 3000° P 

T r = 426° P (all values of thermal conductivity) when 
T e = 5000° P 

The theoretical temperature distribution in the trailing segment 
of the turbine blade was determined by a three-dimensional analysis 
using equation (62) in which the blade cross section was approximated 
by a rectangle. The temperature distribution was found in two planes 
representing the maximum and the minimum temperatures for the z 
axis; the first plane was located at the side of the rectangle and 
the second was on the median plane through the section. The tempera- 
ture was plotted against radial distance from the wheel axis for 
intervals from the trailing edge to the wall of the coolant passage 
(fig . 7 ) . At the wall, the series converges too slowly to be useful, 
therefore the curve was based on a one-dimensional calculation (equa- 
tion (ll)) and the behavior of the other curves of the family. Sim- 
ilar distributions were calculated and plotted using thermal conduc- 
tivities 120 and 210 Btu per hour per foot por °P (fig. 8) . 

The next step in the analysis was to determine the effect of using 
a rectangular cross section as an approximation of the true blade 
shape. The complexity of a solution using the true blado cross 
section eliminates an exact comparison. A wedge-shaped cross section 
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represents the closest approximation with planes as boundaries and 
this shape was used to indicate the trend. A three-dimensional 
analysis of a wedge-shaped cross section was beyond the scope of 
this report and quite unnecessary because on the outer three-quarters 
of the blade, the metal temperature does not change with the radius 
and the comparison was made on a one-dimensional basis for tempera- 
tures along the center line of a section through the outer portion 
of the blade. The comparison is presented as a plot of metal tem- 
perature against distance from the wall of the coolant passage for 
an equivalent rectangular and wedge-shaped section (fig. 9). 

The thermal conductivity of the metals need in blade construc- 
tion has considerable influence on idle temperature distribution; the 
high- tensile -strength alloys currently used have a relatively low 
thermal conductivity. The effect of thermal conductivity was 
investigated using a one-dimensional analysis of a section through 
the outer portion of the blade for average thermal conductivities 
ranging from. 15 (steel) to 210 Btu per hour per foot per °F (copper). 
Although the upper range of conductivities is impracticable for 
actual use in blade construction, it does indioate the trend and the 
magnitude of temperature changes with changes in thermal conductivity. 
The comparison is presented as plots of metal temperature against 
distance frcaa the wall of the coolant passage for the rectangular 
shape (fig. 10(a)) and the wedge shape (fig. 10(b)). The effect of 
a change in thermal conductivity on trailing-edge temperature was 
determined comparatively using the same values of conductivity for 
equivalent rectangular and wedge shapes. One station on the section 
center line at the trailing edge was used for each thermal conduc- 
tivity (fig. 11). 


RESULTS AMD DISCUSSION 

The theoretical three-dimensional temperature distribution for 
a turbine blade of rectangular cross section is presented for effec- 
tive gas temperatures of 2000°, 3000°, and 5000° F in figure 7. The 
temperature falls off in increasing increments as the distance to 
the cooling liquid decreases. At a constant distance from the wall 
of the coolant passage, the temperature remains substantially constant 
for approximately three-quarters of the blade length. The tempera- 
ture beyond this point falls off rapidly and reaches the rim temper- 
ature at the base. The average temperature at the trailing edge of 
the blade was reduced 7.14 percent below the gas temperature at 
2000° F and 7.6 percent below the gas temperature at 5000°. The 
temperature differential between the side and the center of the 
blade (z * t/2) was 3,6 percent of 0, negligible at the trailing 
edge, and at the maximum point did not exceed 5 percent of the metal 
temperature . 
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The temperature distribution at a gas temperature of 2000° 3? 
for rectangular blades having thermal conductivities of 120 and 
210 Btu per houx 1 per foot per °F is presented in figure 8. With 
increasing thermal conductivity of the blade metal, the temperature 
tends to decrease at the trailing edge and increase near the coolant 
passage. The temperature at the trailing edge is approximately 
32 percent below the gas temperature at a thermal conductivity of 
120 Btu per hour per foot per °F and 36 percent below the gas tem- 
perature at a thermal conductivity of 210 Btu per hour per foot 
per °F. With increased thermal conductivity, the influence of rim 
cooling is more effective and the radial distribution lines fall 
off more rapidly. The tempex-ature differential between the sides and 
the center of the blade decreases with increased thermal conductivity. 
For the higher conductivities, this differential is too small to be 
shown in the figures, inasmuch as it ranges from 2° to 7° F. 

An investigation of a rectangular cross-section blade indicates 
the temperature distribution in a blade having the cross section 
shown in figure; 2 but errors in the temperature distribution 
will exist. The trend will be indicated, however, if the rectan- 
gular cross section is compared with a closer approximation to 
tbe true shape; namely, a wedge-shaped section. Figure 9 presents 
plots of the one -dimensional temperature distribution along the 
center line at the blade tip for a rectangular and a wedge section 
at a gas temperature of 2000° F and a thermal conductivity of 15 Btu 
per hour per foot per °F. The temperature distribution for the 
wedge shape has a slightly steeper slope than that for the rectan- 
gular section and the temperatures at the trailing edge and the wall 
of the coolant passage are lower. (Part of this lowering is due to 
the additional thiolaiess of the cooling surface.) 

The effect of the thermal conductivity of the blade metal on 
one-dimensional temperature distribution is shown in figures 10 
and 11. Figure 10(a) presents the temperature distribution along 
the center line of a rectangular section through the outer portion 
of the blade for an effective gas temperature of 2000° F and various 
thermal conductivities. At low thermal conductivities, the temper- 
ature falls off sharply from the trailing edge to the coolant 
passage. With increasing thermal conductivity, the temperature 
decreases at the trailing edge and increases at the wall of the 
coolant passage, which in general tends to equalize the temperatures 
across the width of the blade. The redaction of the temperature 
at the trailing edge increases with an increase in the thermal 
conductivity. Figure 10(b) shows the temperature distribution 
for a wedge section for the same conditions imposed on the rectan- 
gular section. The general form of the curve is similar to that 
for the rectangular section; however, the temperatures at the 
trailing edge and at the wall of the coolant passage are lower and 
the drop in temperature for a given increase in conductivity is greater. 
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The variation in temperature with an Increase In thermal con- 
ductivity for rectangular- and vedge-shaped sections is shown in 
figure 11. The point selected for comparison was located on the 
center line of the section at the trailing edge. The curves show 
that the "blade temperature at this point may be considerably lowered 
by using metals having a high thermal conductivity and that for the 
same conductivity, the cooling is substantially greater for the 
more representative wedge-shaped section. This effect is partly due 
to the loss of effective area at the cooling rasBage where the thick- 
ness of the equivalent rectangular section is considerably less than 
that of the wedge section. At high values of thermal conductivity 
a representative temperature distribution can be obtained from the 
rectangular section by assuming the mean thickness t equal to the 
value of t 2 for the comparative wedge. The total cross-sectional 
area then is no longer equivalent. The same care in selecting 
appropriate dimensions extends the useful range of the nondimenBional 
chart (fig. 12). 

T e - T 

The nondimensional cooling — at the trailing edge for 

“ *! . 

the rectangle is shown as a function of ab' and qo/ak in fig- 
ure 12. Equation (ll) was used with cosh ae «1. It is seen from 
figure 12 that for ab' = 3, the cooling is less than 10 percent 
for the entire range of q 0 /ak. In order to obtain 50 percent 
cooling even for quite large values of q D /ak, ab * must be less 
than 1. Points for the extreme values of k, 15 and 210 Btu per 
hour per foot per °F, are shown on the figure. Theoretically ab' 
could be reduced by decreasing q^ and b ' or by Increasing k 
and t. Under current conditions, not much can be done about q^ 
and k. The curves then indicate short thick trailing edges Instead 
of long thin ones for good cooling. From the curves of figure 12 
It can be deduced that with materials of low thermal conductivity 
an infinite increase In the coolant-passage heat-transfer coeffi- 
cient relative to tha t obtained with water will have small effect 
on the trailing- edge temperature. 


CONCLUSIONS 

The following general conclusions can be drawn from the analysis 
of the cooling characteristics of the trailing part of a liquid- 
cooled turbine blade: 

1. For blade materials having low thermal conductivity, the 
computed three-dimensional temperature distribution of the equiva- 
lent rectangular trailing section indicates that the rim-cooling 
effect is significant for only a short distance from the blade root. 
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the trailing-edge temperature closely approaches the effective gas 
temperature, and the temperature gradient normal to the surface 
exposed to the combustion gas can generally be neglected. 

2. The three-dimensional and one-dimensional solutions of the 
equivalent rectangular trailing portion give almost identical tem- 
perature distribution in the blade section away from the influence 
of rim cooling; the one-dimensional treatment of the wedge or 
rectangular section is therefore sufficient for most investigations. 

3. For blade materials of very low thermal conductivity, the 
one -dimensional temperature distribution computed with the equiva- 
lent rectangular section is only slightly different from that of 
the more representative wedge section. 

4. The equations for the one-dimensional rectangular solution 
permit construction of a simple chart, which gives the value of the 
trailing-edge nondimensional cooling parameter with any possible 
combination of the design and operating variables considered in the 
analysis , 

5. For the blade shapes considered in this analysis, and with 
materials of low thermal conductivity, an infinite increase in the 
coolant-passage heat-transfer coefficient relative to that obtained 
with water will have small effect on the trailing-edge temperature. 

6. For the blade shapes considered in this analysis, trailing- 
edge temperature may be considerably reduced with a material of high 
thermal conductivity and to a lesser extent by locating the cooling 
passage closer to the trailing edge. 

7. Large reduction in trailing-edge temperature results when 
the blade shape is altered to provide a short thick trailing section. 


Aircraft Engine Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 



APPENDIX A 


CALCULATION OF THE INTEGRATION CONSTANTS 


Tho development of equation (54) leads to a deteimlnatlon of N. 


CL Co 

T e - T r = Tjn ]Tm (A^ cos N^I' cosh tWr* cos 4 Z +San cosh Bjnn l* cos X^y' cos P m a) 
1 1 

The first term will vanish if ; 


where 


cog NqI 1 = 0 


- ( n “ |)rt 


n — Z f 4, ■ 1 1 


/ • 1 \ 

\ n - 2) 


For these values of Jr equation (40l) assumes the form 
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xi=kp w 

Tq - T r - >: Bma 0001:1 %n*’ 000 W - ' 000 Pm z 

n=*l m^l 

Equation. (44) is a Fourier development along y 1 and z. 

The development of equation (35) leads to a determination of X. 


CP CSV 

Z, m &■ *(Awx 000 %1 sinii Unnb' cos P^z-B^^coah En^x^Xn sinXnb' cos P m z) 

* (X 

r O <T 


2 n 9o(^e “ Tj-Ajua oos %x* cosh Numb 1 oos Pmz-Bmn cosh K^nx' cos-Xnb* cos P m z) 
11 


In order to simplify aquation (45) at the boundary condition y' ah* and to solve for 
the constant the second surmnatlon in the left member is equated to the last summation in 
the right member term "by tern 


op co 
*r\ <r^ 

n \ t 


9o ( T e ' T l) V (**W Binl1 00831 1 W t ‘) 008 p m z 000 (46) 

1 1 


QC^ QP^ OQ OP 

<L. m Z, n KBjun^oosh I^mpc'jXn 0111 000 ^m 2 <lo Bmn 00011 Emit*' cosXnh 1 oos Pm 2 (47) 

11 11 

The sums are equated term by term. Therefore, X n is determined by any solution of 


tan Xnb 1 
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The development of equation (36) leads to a determination of P. 


CO 

O 


SL, f / 

2^n -kAjnn cos BnX 1 cosh M^ny’ (Pm Bin Pm. gj -^Bnra cosh ’B mn x l coe Xny 1 (^Pm Bin 

00 OCI f t\ 

= ^Tm ^n -<ii Ajnn cos %i T cosh. M^y 1 (cos P m gj -qi Bum oosh E^nX 1 cos X^y' 


1 1 



or 


tan Pm | 



9i* 

2k 

P - 
% 2 


(50) 


His possibility of determining values of Bmn and A^ .to satisfy equations (44) and 
(46) has been established in textbooks on Fourier and other harmonic series (reference 6, 
pp. 118-121), and it is only necessary that tho values determined define a convergent series. 


The values of A rm are deter min ed by integrating equation (46) between the limits of 
x' = 0 to if a I 1 and z = 0 to , z = |> and substituting tho values previously determined 

for Bq and P m in equations (43) and (50). The integration is accomplished in two steps 
using the function cos N B xdx and cos P u zdz as multipliers. 


For abbreviation, 


i|f = (kHjari sinh Mbnb 1 + q Q oosh Mmnb 1 ) 


(51) 


* * 
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Then 

^ Ann 003 Pmz 
because all the tonus 


x’=0 


oos 2 NnX' ax' = (Jo (T e - Tj) 


rx'=z' 

f 

x'=0 


cos N n x’ dx* 


008 H^X* COS HgX' dx' 

Jo 

vanish vhen a / n. Whan oquatlon (52) is Integrated, and tho limits substituted 

003 L ... . sin ^nx’'!* 1 j"(T 0 - Ti) " X ^ 


V Arjrri 


3B ^ Nx , . sin2Ifnx>'] |(T a -Ti) 1 

2K n KnX + 2 3 <4 — K Bln 

*<x»=0 ** ^ x 1 =0 


or 


V 4m oos P m z (V + » 2 ^ — ~ ~ sin N n i' 


2Hn / 

Whotn, tho value % as determined in equation (43) is substituted 

n-1 


(52) 


(53) 


(54) 


fJU °ob Pni . i (T 0 . t 2 ) hi)?! „ Sa 
When oquation (55) is integrated in teimo of s using tho multiplier cos P u zdz, 


(55) 


r Z' 


, i T . 

V -Amn. oob 41 P m zdz & p n ooe p^ziz 

Jz=0 ^z=0 


(56) 


co 

H 
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"because sll the terns 


i 2 


| cos P m z cos P u zdz 
,,'o 

vanish when u /n. The integration of equation (56) gives 


Substituting limits, 


.. t 

-2=o 


Vi sin 2Pa z 

V 5^=- P W Z + 


2P m | x »‘ 


Pn 

Pi 


z=^0 


c h h 


t 

-2 


V Anm. (?r 


sin 2P m 2pn sin P n | 
v 2 + 2P m / P a 


(57) 


(58) 


ro 

CO 


or 


Pn sin Pm g- ^ 


(1 + 

V 


2P m | 


In like mann er Bum is determined from equation (44) and tho equation for is 

1 


Bmn = 2£n 


sin Pn !■ 


P n | /' sin 2P E I' 

cosli B wn I 1 ( 1 + r~ 

2P n 2 


(59) 


( 60 ) 
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■whoro 




2 


( 3 ?e “ *r) 


sin T^b’ 1 

Anb 1 f n sin 2Vo'\ 

V 1 + 2^* ) 


( 61 ) 
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APPENDIX B 


EVALUATION OP THE TEBM3 WHERE m > 1 

The temperature distribution in the turbine blade under consid- 
eration is determined neglecting the smaller terms when m > 1 and 
using equation (37) in the form 

OS 

9 = cos P]_z n (Ain cos N n x’ cosh M-) n y 1 +Bi n cosh Ri n x' cos A n y 1 ) 

1 (62) 


This equation will be valid for all practical purposes when the 
thickness of the blade is small compared with the width and the 
length. It is necessary, in demonstrating this statement, to expand 
the double series given by equation (37) and inspect the terms con- 
taining the function cos P m z. The expansion of equation (37) with 
the subscripts in the order mn assumes the form 

9 = (An cos l 1 Mjjy* cos P^z+Ag]_ cosh Mo^y* cos Pgz + ...) cos N^x* 

+(Ai 2 cosh cos P-^z + Agg cosh Mggy* cos P 2 z + ••») cos Ng 1 ’ 

+. .. + (Ajj^ cosh M^y 1 cos Pqz+A^ cosh Mg^y* cos P 2 z + . ..) cos N^x* 

+(Bn cosh E^qx* cos Vjz. +B 2 ^ cosh. cos P 2 z + .. *) cos A]y* 

+(B^2 cosh B^gx' cos P^z+Bgg cosh Bggx* cos Pgz + •••) cos Agy 1 

+ . . . + (Bpn cosh Ei-nX 1 COS P;j_Z cosh Bg^X 1 cos P n z + . . . ) cob A^y 1 

( 63 ) 

If all the other terms in the series are arranged so they become 
coefficients of the function 00 s P m z, 
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Q = An cosh Mu? 1 (°ob P x z + ^ ^ 008 P 2 Z + --j 088 K l*‘ 


/ A 22 COSh Mgoy 1 \ 

+ A 12 cosh %27' |cos P X z + ^ M ^ yt cos P2Z + . • y cos N 2 x* + . . 


a . „ A2n cosh M?ny 1 _ \ 

+ A ln cosh ^y* ^cos P x z + CQsh M ^ yt cos P 2 z + . . . | cos %x» 


i 


_ . _ , / _ B 2 i cosh E2I3C’ \ 

+ B 1l cosh H u i* ^cos P x z + Bn " cosh Bjjx t 008 P 2 Z + • • -J 008 *a7* 

( Bpp cosh Rppx* \ 

R 1 2 x ’ <° os P l z + B^ '~ c8sh R^x» 888 P 2 Z + • • •) 888 *2^ 


+ B X2 cosh 


/ Bgn cosh Rj>nX 1 \ 

+ B^ cosh R^x 1 ^cos P x z + B ^ ' obsh R-^x* 008 P 2 Z + • • • i 008 

(64) 


From equation (64) It can ho seen that if the coefficient of the 
function cos P 2 z in any series Is negligible compared with unity, 
the terms subsequent to cos P x z may be disregarded for practical 
application of the equation. Because the series is convergent for 
ascending values of n, it Is only necessary to inspect the coef- 
ficient of cos P 2 z for a value of n = 1. 


When the coefficient of cos P 2 z in the series Is equated 

to 0 and the values of ^ determined by equation (59) sire substi- 
tuted. 


sin P 2 \ f sin 2P X 

t " ^ ^ 

p 2 2 " \ ^1 2 " / (Mil i sinh M xx b’ +qo cosh Mn b 1 ) 

sin Px | 7 sin 2P 2 |\ 8inil M 21 V cosh M 21 b*) 


cosh Mgjy* 
cosh Mi i y* 
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or 


ein P 2 \ ^ sin 2P 1 f\ 


* 2 I \ 


( 1+ 


2P 1 2 J (kMll tanh MiiV + (jo) cosh Mi i b 


sin Pi % f sin 2 P 2 g 
r — [1 + — 


cosh Msir’ 

^ (kM21 tanh M 2 ib J + g^) cosh M 21 V cosh Mny T 


’ 1 2 


\ 


2F2 I 


1 


( 66 ) 


For the blade under consideration, 
V = 0.055 ft 
k = 15 Btu/ (hr) (ft ) (°F) 

1* = 0.3333 ft 

a± = 222 Btu/ (hr) (esq. ft)(°F) 

= 2370 Btu/ (hr)£oq ft)(°F) 
t = 0.01 ft 


From equation (50) 

P L = 53.74 

P 2 = 632.99 

From equation (43) ...... - 

I*! = 4.64 ft" 1 
From equation (38) 

Mu = 53.94 ft" 1 

M 2i = 632.99 ft" 1 

33ie following circular and hyperbolic functions apply for these 
values : 
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Bin P x | = 0.2654 

sin Po £ = .0234 

^ 2 

sin 2P 1 | = 0.5119 
Bin 2P 2 | = .0467 

tanh M i;l V = 0.995 
tanh » = 1 


cosh M-j^h ' = 
cosh Mg-^h* = 
cosh Mny' = 
cosh Mg^y* = 
cosh M^y' = 
cosh Mg^y' = 


9.7425 

6.5907 X 10 14 
1 when y* = 
1 when y ' = 
7.4513 when 
2.7822 x 10 13 


0 

0 

y* = .05 
when y' = .05 


If the given values in equation (66) are substituted,, the following 
coefficients are obtained: 

y 1 = 0 0 = -5.738 X 10 ~ 17 

7 ’ = 0.05 0 = -2.142 x 10" 4 

y* = b» 0 = -3.882 x 10 “ 3 

The values of f> obtained for the blade under consideration show 
that the coefficient of the cos Pgz term has a maximum value at 
y 1 = h 1 and is then equal to 0,4 percent of the first term and as 
such may he disregarded for practical applications . 


The effect of the second term of the series can he inves- 
tigated hy equating the coefficient of the cos Pgz term to x and 
the previous process repeated. 


sin 


p 2 I /f sia lr\ 

”t l 1 + 7T t 1 


I t 1 OT , t f 

P 2 2 \ 2Pi g* j cosh E X1 I ’ cosh E 21 x* 

~ sin P X | / sin 2Po £\ °° 8h *21* ’ 


1 + 


4 rrrt t\ COSh EotI ' 

sin 2P 2 2 \ 


(67) 


p iJ V 


2P-; 


* J 


The application of equation (67) to the blade under consideration 
requires the following values: 


From equation (48) 


\ L = 25.64 ft" 1 
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Prom equation (39) 


r?ll = 59.54 ft" 1 

S 2 l = 633.51 ft" 1 


With the values of B^ and 7' used, the cosh I Wi Z ' and 

g BTTTTl Z * , . 

cosh B^x 1 may he approximated hy g and equation (67) 

becomes 


sin Pg Z sin 2Px ^ 


r * 


P„ % 
2 2 


2P 


sin Px \ / sin 2P 2 £\ 

r i| ( + 2p a I ) 


1\ 

I 1 ) 

LL^. e (Z'-x')(B 11 -B 21 ) 


( 68 ) 


When the given values in equation (69) are substituted for a 
range of values of x 1 , the following values of T are found: 


X 1 

(ft) 

I Distance from rim 
(ft) 

r 

0.005. 

0.3335 

A 

.3283 

.01 

-0.4643 X 1C" 4 

.3333 

.005 

- .8227 X 10" 3 

.3383 

0 

- .1451 X 10" 1 


The coefficient of the cos P m z term increases as the rim is 

approached and reached a value of 1,5 percent of the first term at 
thiB point. For purposes of this report, the temperature will not 
be calculated within 0.005 foot of the rim and the cooling liquid. 
With this limitation, the values of the second and subsequent terms 
are negligible . 
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Figure 6. - Correction for heat received by trailing edge 
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Blade length, in. 

(b) Effective gas temperature, T e> 3000° F. 

Figure 7. - Continued. Three-dimensional temperature distribution in rear part of tur- 
bine blade at thermal conductivity k ot 15 Btu/(hr)(ft)(°F). 
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Fig. 7 c 



Figure 7. - Concluded. Three-dimensional temperature distribution in rear part of tur- 
bine blade at thermal conductivity k of 15 Btu/(hr)(ft)(°FJ. 
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Figure 8 . - Three-dimensional temperature dlsti Ibution In rear part of turbine blade at ef- 
fective gas temperature T e of 2000 s f; 
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Fig. 8b 



Blade length, in. 

(b) Thermal conductivity, K, 210 Btu/( hr ) (ft) (°F) . 


Figure 8. - Concluded. Three-dimensional temperature distribution in rear part of turbine 
blade at effective gas temperature T e of 2000 w r . 
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Figure 9. - One-dimensional analysis of the effect of shape on prevelanet blade tem- 
perature. Effective gas temperature, T e , 2000° F; thermal conductivity k, 15 Btu/ 
(hr) {ft) (°F) . 
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(b) Wedge-shaped section. 

Figure 10. - Concluded. One-dimensional solution for prevalent blade temperature at 

effective gas temperature T e of 2000° F. 
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Fig. 12 
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Figure 12. - Nondliwneional cooling at the trailing edge T«-T/T«- T. as a function of d i mention I eee param- 
eter* ab' and q 0 /aK when a z V*q]/Kt . 
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